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Abstract 

The order of the Coulomb-Higgs transition in the U(l)-Higgs model with 
unfrozen modulus of the scalar field is studied. Large lattices (up to 24^ in 
one case) and high statistics are used. We fix /? = 1.15 and explore specially 
a region of A- values where metastability is observed. We study the thermody- 
namical limit of several observables, in particular, the latent heat, the specific 
heat, the decrement of the free energy between the maxima and the central 
minimum of the two-peaked histogram, the Binder cumulant and the displace- 
ment of the critical coupling with the lattice size. The results point towards a 
second order transition for A > 0.005, while for smaller values of A the strong 
metastability growing with the lattice size seems to derive from a first order 
character. 
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1 Introduction 



The Higgs mechanism is an essential part of the present day formulation of the 
Standard Model. The U(l) gauge model coupled to scalars is a simplification of the 
gauge-Higgs sector of the SM, which keeps the unsolved problem of defining non 
perturbatively a non-asymptotically free field theory. 

The U(l)-Higgs model has been widely studied previously. One of the main 
objectives has been to determine the order of the Coulomb-Higgs transition, both 
in the cases of frozen and unfrozen modulus of the scalar field, by using analytical 
(mean field, one-loop effective potential, etc) and numerical techniques (see || 
and references therein). 

From those analyses the situation could be defined in the following way. For 
large values of A the transition seems to be second order, while for small values 
it looks first order. These results are based on Monte Carlo studies of the system 
searching for metastabihties or single/double peaked histograms. Due to computa- 
tional limitations, these calculations had been made with small lattices and short 
statistics. Here we carry out a study with much larger lattices and high statistics 
in order to approach the thermodynamical limit in a more reliable way, obtaining 
results qualitatively consistent with the previous ones. 

However, in those works the conclusion that the transition is first order has been 
obtained by considering the presence of a double peak for a given volume V (or 
observing metastability). As we will show this is not correct because even in this 
case, when V is increased, both peaks approach, and the latent heat disappears in 
the thermodynamical limit, obtaining in this way a second order transition for A 
values much smaller than previously considered. 

2 The model and the Coulomb-Higgs transition 

The three parameter U(l)-Higgs model is described by the action 

r,^<p r,fi 



In the A ^ oo limit, |$| ^ 1 and the action simplifies to 

S = -f3Y^ 3fJf/,,^, -kY, 3?$.f/.,^$.+^ (2) 



r,fj.<u r,fj. 



The phase diagram of that restricted version was considered first in P] and has 
been discussed by us in and 0. We discuss here the global aspects of the phase 
diagram in the fixed modulus case (See Figure |l|). Point A is the pure compact 
U{1) phase transition, a well established first order point even though 
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Figure 1: Phase diagram of the model with fixed modulus. 

this transition for small lattices seems to be second order. We point out that some 
authors have cast a new hue on its nature, pointing out the possibility of this 
transition to be second order. 

As we showed in EI, particularly in the neighbourhood of the triple point C, the 



line joining A with the latter is also first order. The line CD was studied by us [|10 
and is a first order line ending in a second order point D. The precise location of D is 
K = 0.5260(9) and f3 = 0.8485(8), with measured critical exponents compatible with 
the classical (Mean Field) values a = 0, (3 = 1/2,7 = 1, = 1/2. The line BC is 
more controversial. The difficulty to identify the order of the Coulomb-Higgs phase 
transition was made apparent because of the large finite size effects. The (3 ^ oo end 
of that line is the X-Y model phase transition, a well established second order one 
with Mean Field exponents, which has induced researchers to think the line to be a 
prolongation of its end point, and as such, second order |ll|. Yet, the relationship 



of the model to the q > 1 version and to the model reported in ||T2|, and [g 
points towards a possible first order transition. However that relationship is based 
on perturbative arguments (expansion around f3 = oo) and might not be applicable. 

The difficulty to study directly the A — oo limit has lead us to study the problem 
at finite, variable A. This had been done earlier 0, and we intend to improve on 
the statistics and the lattice sizes. 
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3 Simulation and observables 



We have fixed [3 = 1.15, which is on the Coulomb-Higgs side of the first order 
Confining-Coulomb phase transition. 

If we consider larger values of /3, the system has a larger correlation length, and 
then in order to approach the thermodynamical limit, L ^ ^, we need very large 
volumes. Also, a tricritical point along the Coulomb-Higgs transition, where the 
transition order changes, is not expected. 

Then, we select some value of A and search for the "critical" k, i.e. we are 
looking for the Coulomb-Higgs transition. As is known, for small enough values 
of A the simulations on this transition show clear signals of met ast ability, giving 
place to two-peaked histograms, while for large values the metastability disappears. 
Of course, the simulations are made on finite lattices and then a conclusion on 
the order of the transition cannot be extracted directly from that behaviour. We 
have selected an intermediate region of A values, ranging from the point where the 
two peak signal practically disappears for our lattices sizes (A = 0.3) to the point 
where the metastability is so strong that it makes the work with our computational 
disponibilities difficult (A = 0.003). The total set of A- values that we have used is 
0.003,0.005,0.01,0.03,0.1,0.3 on lattices 6^ 8^ 12^ 16^ and 24^ (depending on A), 
with statistics ranging between 10^ and 10^ Monte Carlo iterations per measurement, 
on workstations, on parallel machines at the Edinburgh Parallel Computing Centre 
and on our custom, 64 T800 processor computer RTN |T^. We have used an over- 



relaxed Metropolis update method, with a multicanonical procedure to accelerate 



the fiip-fiop rate in some case. We use the spectral density technique in order 
to improve the determination of the transition point, and the jack-knife method in 
the error estimation. 

We have observed the different energies one can define from the action. In 
particular we have paid special attention to the link energy -Eunk = (5>t/$) as we 
have fixed the value of f3 and the transition is almost parallel to the (3 axis. All the 
observables we will define hereafter will be referred to this energy. 

In the range of A and L studied, the histograms of -Eunk present the following 
features (examples can be seen in fig. ||): 

• Two-peak structure. 

• Asymmetry of the peaks, with a narrow low-energy peak and a broad high- 
energy one. 

• Strong size dependence, consisting on a narrowing of the gap and width of the 
peaks as the lattice size L increases, caused by a much faster displacement 
towards lower energy of the broad peak. 

It is known that, due to finite size effects, two-peaked histograms on finite lattices 
can become single peaked in the thermodynamical limit and vice versa. For instance, 
the transition in the q = 4 Potts model is known analytically to be second order 
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but the simulation on finite lattices gives a double peaked structure [0. On the 
contrary, single peaked histograms on small lattices can become double peaked when 



the lattice size increases, as in IK 



Then an accurate study of finite size effects has to be made. We will use several 
independent Finite Size Scaling methods, which allow to determine correctly the 
order of the transition whenever one uses large enough lattices (but not so large as 
if one only looked for metastability signals). Now we will introduce the different 
quantities which we have considered in order to extract our conclusions. 



First, coming back to the histograms, in ref. [0 an interesting method is pro- 
posed. They explain how in a two-peaked structure the increment between the 
minima and the local maximum of the free energy is directly related to the analo- 
gous increment for the logarithm of the histogram of the energy. The idea is that if 
the transition is first order the gap has to be more pronounced for larger L, while if 
it is second order the depth of the gap does not grow with L. Our procedure will be 
the following. From the histogram of -Eunk obtained in the simulation, we calculate 
a new histogram where the height of the two peaks is the same by using basically 



the spectral density method [Q. Then we take the logarithm of it and change its 
sign. In this way we obtain a figure analogous to the free energy (see fig. Here 
we can measure directly the depth of the gap, let us call it AF. A AF growing with 
L indicates a first order transition, while if it is constant we have a second order 
one. 
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Second, we have measured the specific heat 



Its maximal value should scale with an exponent a/i>. It will be d for a first order 
transition and for a mean field second order one (in this case c™"^ will be a constant 
plus logarithmic corrections). 

Third, we use the maximum of Cy as definition of critical kappa Kc{L,\). The 
scaling of ^^{L) with L allows us to extract v. It will be z/ = l/rf if the transition is 
first order, and t/ = 1/2 if we have second order with mean field exponents. 

Fourth, we consider the latent heat at Hc{Li\)-i AE", as the difference between 
the positions of the maxima for a fit to each peak separately to a cubic spline, after 
the histogram has been shifted to the apparent critical point by the spectral density 
method. A cubic spline has been preferred to other functional shapes because of its 
ability to reproduce a maximum and accommodate the mixed states, whose influence 
on the histogram is otherwise difficult to account for. If we have AE 7^ in the 
y — > 00 limit we have a first order phase transition. 

Fifth, we have computed the Binder cumulant 

B^l- (4) 

If the energy distribution is single peaked (and then second order) in the thermo- 
dynamical limit, this quantity can be computed exactly and becomes 2/3. 



4 Results 

First of all, let us comment on which have been our limits on A and L. 

For the biggest value of A studied, A = 0.3, it was impossible to distinguish the 

two peaks on the histogram, and then we could not make the type of treatment 
described above. For A = 0.1 two peaks arc visible but the gap is so small that it is 
difficult to measure with precision observables such as AF and A.E; for this value 
of A we have only simulated L = 6, 8. 

On the other hand, for the smallest value A = 0.003 the metastability is so 
strong that we have only been able to use L = 6, 8, 12. Besides, the simulation in 
the L = 12 have been made fixing a value k = 0.2703 and making two different runs, 
one starting from a cold configuration and the other one starting from a hot one. 
In this case (L = 12) no flip was observed; then we can only say that the two peaks 
are rather separated, but a precise estimation of the latent heat and of the other 
observables is not possible. In particular, the latent heat is greatly affected by a 
displacement in k, (because the broad peak is much more sensitive than the narrow 
one, as the displacement in the peak will be proportional to the squared width) and 
a bad estimation of Kc would give place to an erroneous measurement. 
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Figure 3: Latent heat as a function of at different values of A. From top to 
bottom, A = 0.003, 0.005, 0.01, 0.03, 0.1. 



A 


L = 6 


L = 8 


L = 12 


L = 16 


L = 24 


0.003 


7710(490) 


9810(1150) 








0.005 


1820(80) 


2400(70) 


3850(170) 


4260(890) 


7760(810) 


0.01 


577(57) 


652(53) 


660(28) 


897(105) 




0.03 


185(41) 


189(17) 


172(15) 


186(12) 




0.1 


63(2) 


58(3) 









Table 1: Specific heat for the different values of L and A. 



For the intermediate values of A things are more feasible and we have used 
L = 6, 8, 12, 16 for A = 0.005, 0.01, 0.03 and also L = 24 for A = 0.005. 

Let us now present the results for the different observables. 

The study of the latent heat suggests an extrapolation to zero in the thermody- 
namical limit for A > 0.005 as shown in fig. ^. 

The specific heat cy is constant with L for the largest A-values while it grows 
for the smallest ones, as shown in table |l]. In the presence of a stable double peaked 
structure (i.e. such that the peaks do not approach when changing L), cy scales 
with the maximum possible speed, i.e. like V , indicating a first order transition. 
However in our case the peaks are not stable and then the scaling will be different. 
For a trivial second order transition in d = 4, it scales only logarithmically (a = 0), 
then it must be almost constant. The results indicate a crossover between this two 
behaviors at A ~ 0.005. 

For the values of A where we have more measurements, i.e. A = 0.005, 0.01, 0.03 
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Figure 4: Binder cumulant at A = 0.005 as a function of L^^'^. The circle corre- 
sponds to the value 2/3. 

we have fitted Kc{L) as a function of L: 

K^{L) = K,{oo) + AL-^I\ (5) 

The three parameter fit gives v = 0.40(1), 0.43(1), 0.44(1) respectively at those 
values of A. These values of u are closer to the mean field one, u = 1/2, than to the 
expected in a first order transition, u = 1/4. However, those apparent zz-exponents 
would probably derive to one of the two values (1/2 or 1/4) in the thermodynamical 
limit. 

The behaviour of the Binder cumulant has not given any conclusive result. In 
principle, it should approach the value 2/3 with corrections of the order of L'^^^^'^~'^^ 
for a second order transition or a value different from 2/3 with corrections as L~'^ 
for a first order one. In practice, we obtain some value close to 2/3 but neither the 
precision is enough nor the scaling is as cited. In fig. ^ we show the results for 
A = 0.005. We have a good scaling but with an exponent —2.5(2) instead of —4, 
which would be expected both in the first order case and in the gaussian second 
order one. The fit gives a value 0.64(1) in the infinite volume limit. However this 
difference from 2/3 is not significant as the error is large and the scaling is still not 
as expected in the thermodynamical limit. 

Finally, in table |^ are shown the increments of the free energy AF. Again we can 
see that the behavior changes at A ~ 0.005; for larger values AF is constant with L, 
while for smaller ones it grows. A look at fig. ^ shows the different kind of behavior 
for A = 0.01 and A = 0.003; while in the first case the gap tends to disappear when 
increasing the lattice size, in the second one it tends to become deeper. 
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A 


L = 6 


L = 8 


L = \2 


L = 16 


L = 2A 


0.003 


4.0(4) 


5.0(7) 


large 






0.005 


2.3(5) 


2.1(4) 


3.1(4) 


3.4(4) 


3.5(5) 


0.01 


1.7(4) 


1.9(4) 


1.3(4) 


1.9(4) 




0.03 


0.6(2) 


0.8(3) 


0.6(2) 


0.7(2) 





Table 2: AF for the different values of L and A. 



5 Conclusions 

We have investigated the order of the Coulomb- Higgs transition at fixed [3 = 1.15 
and for different values of A. The analysis of several observables indicates a crossover 
at A ~ 0.005, the transition being first order for smaller and second order for larger 
values of A. We remark that in order to extract this conclusion we have needed 
to consider different lattice sizes and, even though most of them present a two- 
peak structure, the latent heat goes to zero in the thermodynamical limit (at least 
for A > 0.005). Previous works had estimated this crossover at A ~ 0.3 [Tl|] or 
A ~ 0.02 |I| for (3 = 1.5. However for larger (3 a weaker character is expected 
(the correlation length grows and we approach the second order transition of the (j)'^ 
theory). Then for (3 = 1.5 the crossover should occur at a smaller A- value than for 
f] = 1.15. Therefore we have lowered the estimation of the values of A for which the 
Coulomb-Higgs transition changes order. 

Acknowledgments 

We are very grateful to the RTN collaboration for the use of their transputer machine 
where we have run most of our simulations. D.I. acknowledges the EPCC and the 
Deparment of Physics of Edinburgh, as well as the Ministerio de Educacion y Cultura 
(Spain) for a fellowship. This work has been partially supported by CICYT under 
contract numbers AEN96-1634 and AEN96-1670 and by the EU under contract 
number ERB-CHRX CT 92-0051. 

References 

[1] D. Espriu, J.F. Wheater, Nucl. Phys. B258 (1985) 101. 

[2] H.G. Evertz, K. Jansen, J. Jersak, C. B. Lang and T. Neuhaus,Nucl. Phys. B285 
[FS19] (1987) 590. 

[3] E. Fradkin and S.H. Shenker, Phys. Rev. D19 (1979) 3682. 

[4] A. Cruz, A. Tarancon, L. A. Fernandez, A. Munoz Sudupe, and J.J. Ruiz- 
Lorenzo, Nucl. Phys. B (P.S.) 34 (1994) 638. 



9 



[5] A. Cruz, A. Tarancon, L. A. Fernandez, A. Munoz Sudupe, and J.J. Ruiz- 
Lorenzo, Int. Jour, of Mod. Phys. C, Vol. 5 No. 2 (1994) 343. 

[6] M. Creutz, L. Jacobs and C. Rebbi, Phys. Rev. D20 (1979)1915. 

[7] B. Lautrup and M. Nauenberg, Phys. Lett. B95(1980)63. 

[8] L. A. Fernandez, A. Munoz Sudupe, R Petronzio, and A. Tarancon, Phys. Lett. 
B267 (1991) 100. 

[9] J. Jersak, C.B. Lang and T. Neuhaus, Phys. Rev. Lett. 77 (1996) 1933. 
[10] The RTN Collaboration, J.L. Alonso et al. Nucl. Phys. B405 (1993) 574. 
[11] K. Jansen et al, Nucl. Phys. B265 [FS15] (1986) 187. 

[12] L. A. Fernandez, A. Muiioz Sudupe, J.J. Ruiz-Lorenzo and A. Tarancon, Phys. 
Lett. B312 (1993) 305. 

[13] RTN Collaboration, Computing in High Energy Physics '92 CERN 92-07, 
Geneva (1992) 353. 

[14] M. Falcioni, E. Marinari, M.L. Paciello, G. Parisi and B. Taglienti , Phys. Lett. 
B108 (1982) 331; A. M. Ferrenberg and R: Swendsen, Phys. Rev. Lett. 61 1988) 
2635. 

[15] M. Fukugita, H. Mino, M. Okawa and A. Ukawa, J. Phys. A. Math. Gen. 23 
(1990) L561. 

[16] J.L. Alonso, J.M. Carmona, J. Clcmcntc Gallardo, L.A. Fernandez, D. Iniguez, 
A. Tarancon and C.L. UUod, Phys. Lett. B 376 (1996) 148. 

[17] J. Lee and J.M. Kosterhtz, Phys. Rev. Lett. 65 (1990) 137. 



10 



